A chiral lagrangian with broken scale: 
a Fortran code with the thermal contributions of the dilaton field 
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The Chiral Dilaton Model is a chiral lagrangian in which the breaking of scale invariance is 
regulated by the expectation value of a scalar field, called dilaton. Here we provide a Fortran 
code P, 0], as a tool to make calculations within the Chiral Dilaton Model at finite density and 
temperature. The calculations are improved respect to previous works 0,131 by including the thermal 
contributions of the dilaton field. 



PACS numbers: 21.65.Mn, 12.39.Fe, 21.65.Cd, ll.lO.Wx 
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I. INTRODUCTION 

In the next years several Heavy Ion Collision (HICs) 
experiments at energies of the order of a few ten A GeV 
(as e.g. the ones proposed at facility FAIR at GSI [5|, 
at RHIC (Brookhaven) and at the Nuclatron in Dubna) 
will probably start their activity. In these experiments 
the Equation of State (EOS) of matter will be tested at 
large density and/or temperature. It is therefore very 
important to provide, through theoretical investigations, 
a map of the "new" regions which will likely be explored 
by the new experiments. 

The Chiral Dilaton Model (CDM), developed W the 
group of the University of Minnesota jj, ^, d, Q and 
largely discussed in Refs. [1, 01, is a chiral hadronic model 
where chiral fields are present together with a dilaton 
field which reproduces, at a mean field level, the break- 
ing of scale invariance taking place in QCD. The presence 
of two not-vanishing condensates, the chiral and the dila- 
ton condensates (the latter is connected to the gluon con- 
densate), allows to study the restoration of both chiral 
symmetry and scale invariance and the interplay between 
them. 

In this work we provide a Fortran code P, 0] as a 
tool to explore the properties of the CDM. In partic- 
ular the code allows to compute the EOS, the masses 
and the mean values of the fields in a wide region of the 
density-temperature plane (O.Olpo < Ps < S.Olpo with 
po 0.15/m~3 and < T < 245 MeV) and in a wide 
range of isospin asymmetries (0.3 < Z/A < 0.5). The 
code allows also to compute neutron star matter at T=0 
(with electrons and muons in /3-equilibrium) and T=0 
pure neutron matter. A very small execution time is en- 
sured since the code makes a cubic spline interpolation 
of data previously computed. 

There are two important improvements present in the 
code respect to the calculations provided in Refs. [1, 
The most important one is that the thermal fluctuations 
of the dilaton field are taken into account. This allows to 
study more in detail the regions of density and temper- 
ature where scale invariance is restored. The second im- 
provement concerns the technique used to compute the 
thermal averages, which is more reliable than the one 
used in Refs. [3, [j] because it does not resort to any ap- 



proximation. 

The paper is structured as follows: in Sec II the la- 
grangian of the CDM is described, in Sec III the tech- 
nique used to compute the thermal averages is presented, 
in Sec. IV the procedure used to make calculations of the 
model is described. Finally in Sec. V we conclude the 
paper by showing some applications of the code. 



II. THE CHIRAL DILATON MODEL 



The lagrangian of the CDM reads: 
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where the potential is: 
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Here a and tt are the chiral fields, 4> the dilaton field, 
the vector meson field and the vector-isovector meson 
field, introduced in order to study asymmetric nuclear 
matter. The field strength tensors are defined by F^i, = 
dpUJv-diyOjp, B^,, = d^h^-d^h^. In the vacuum = 0o, 
a — and tt = 0. The uj and p vacuum masses are 
generated by their couplings with the dilaton field so that 



^y^f 00 and TOp 



G^J^iia- Moreover C = ^o/o-q, B 
and 5 are constants and e'l is a term that breaks explicitly 
the chiral invariance of the lagrangian. Finally, the values 
of the parameters used in calculations (listed in table 1 
of Ref. ^) were determined in Ref. [7] by fitting the 
properties of nuclear matter and finite nuclei. 
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III. THERMAL AVERAGES 



In finite temperature mean field approximation, when 
the temperature is large enough to become comparable 
to the masses of the meson fields it is mandatory to take 
into account the thermal fluctuations of those fields. In 
Ref. [3] we considered the thermal fluctuations of the 
chiral fields and of the vector mesons, but the thermal 
fiuctuation of the dilaton field was not considered due 
to the large value of the glueball mass (about 1.6 GeV 
in the vacuum). In the present work we extended the 
previous calculations by including also the fluctuation of 
the dilaton field. 

As already discussed in Ref. [1] , it is not trivial to com- 
pute the thermal averages of a quantity depending on the 
fluctuating flelds, in particular when those fields do not 
appear in polynomial form. Due to the complicate struc- 
ture of the potential ([3]), almost all quantities depending 
on the scalar fields are difficult to average (notice that 
the terms containing only the vector fields are not prob- 
lematic) . 

In Ref. [3| we presented two different techniques to 
compute the thermal averages: the technique developed 
by the authors of the CDM [1, which we mainly 
adopted in Refs. [1 i], and a technique developed in a 
more recent paper by Mocsy and Mishustin [lol ]. 

The first technique jl, [§] was developed to treat in a 
compact way the chiral invariant term (cr^ + n'^), which 
appears frequently into the lagrangian. To this purpose 
two approximations were introduced: the fluctuations of 
the chiral flelds are assumed to be equal and the term 
ctAct, supposed to be small, is expanded up to the fourth 
order. The second technique |l0|, instead, does not resort 
to any approximation and the thermal fluctuations are 
taken to all orders (notice that this request is not satisfied 
in the previous technique, because the expansion of the 
term a Act is truncated to the fourth order). Moreover the 
technique of Ref. [To| is a general procedure extensible to 
a generic number of fields, where each thermal fluctuation 
is treated as a independent quantity (this allows to treat 
separately each isospin component of the pion, enabling, 
in this way, to compute the isospin splitting of the pion 
mass). 

For these reasons and since we are interested to include 
in the calculation the thermal fiuctuation of the dilaton 
field we adopt in the present calculation the technique 
by Mocsy and Mishustin (T^]. Extending the averaging 
procedure of Ref. by including the fiuctuation of the 
dilaton field and treating the isospin components of the 
pion as independent of each other, we obtain that the 
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where 

F(z„ (A^)) = (2^(A2))- V2 exp (-^) (4) 

is a gaussian weighting function depending on the ther- 
mal fiuctuation (Af ) of the i-th scalar field. Notice that 
this averaging procedure ensures the equivalence between 
the resolution of the field equations and the minimization 
of the thermodynamical potential (this equivalence was 
only approximate with the procedure developed by the 
authors of the CDM), as it has been proved in Ref. p^ . 



IV. FINITE TEMPERATURE CALCULATIONS 

In this section we summarize the procedure adopted to 
make calculations at finite temperature. 

The symmetries of the lagrangian lead to the conserva- 
tion of the baryon number and of the isospin. Depending 
on the thermodynamical ensemble used, the conservation 
of these charges can be expressed by fixing their chemical 
potentials or their densities. Since we use the canonical 
ensemble we fix the baryon density ps, the total charge 
per baryon Z /A (that is equivalent to fix the total isospin 
density) and the temperature T. 

Differently from what done in Ref. [3| (where the mean 
values of ct and <j) were evaluated by minimizing the ther- 
modynamical potential), here the calculation consists en- 
tirely on the numerical resolution of a system of equations 
in the following variables: 

- the nucleon chemical potentials /ip and /i„. 

- the thermal fiuctuations (Acj)^), (Act^) , (ttq), (tt^), 
(ttI), (Auj^Aiv^), (A5o^_A6(;), (A6+^A6^), (Afe_^A6'i) 

- the mean field values (f), ct, loq, bo. 

As it will be shown in the following subsections, 2 of the 
15 equations come from the conservation of the baryon 
and isospin charges , 9 are the self-consistency relations 
for the thermal fluctuations and the remaining 4 are the 
field equations. 



A. Conserved charges 

In order to fix the total baryon charge and the total 
electric charge one has to solve the following 2 equations: 
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where rrii is its effective mass, fj,* its effective chemical 
potential and rii the degeneracy factor: 



_ 1 for the 0, a, ttq, 7r+ and 7r_ mesons 
* ^ ' — 3 for the lo^, b^^, and mesons ^ ' 



In order to compute the thermal fluctuation of a meson 
field (fTTjl one needs the value of its effective mass, which 
in turn depends on the thermal fluctuations and on the 
mean values of all the flelds entering the expression ([7]). 
This leads to a system of 9 self-consistent equations which 
reads: 



where the index i runs on the scalar and vector flelds, 
Pi(m,yU,T) are density integrals, and the effective chem- 
ical potentials of the nucleons which enter the thermo- 
dynamical integrals are related to the standard ones as 
follows: 
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The effective masses of the mesons are given by: 
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where the minus sign stays for the scalar mesons and the 
plus sign stays for the vector mesons. Finally the nucleon 
mass reads: 



'AT 



(8) 



Notice that the effective masses depend on the mean 
values and on the fluctuations of the flelds which are 
needed to compute the thermal averages ^ and ([8]). 

The chemical potentials associated to the charged iso- 
vector mesons are related to the chemical potentials of 
the nucleons by the relation Q: 



A^7r+ pp^ Pp pn 
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and the effective chemical potentials, which enter the 
thermodynamical integrals are given by: 

= Pp+ = Pp- Pn- (10) 



B. Self-consistency relations for the thermal 
fluctuations 



The thermal fluctuation for a generic meson field is 
jiven by: 



^ I '^^e-(0„(A|))exp[/3(e-(0,,(A^^»-p*)J-l(13) 



where the indexes i and j run on the scalar fields 0, cr, 
ttq, 7r_|_, 7r„ and on the vector flelds w^, 60^^, b+^, b-^. In 
the case of isospin symmetric nuclear matter, the thermal 
fluctuations of the isospin components of the pion and of 
the p— meson are equal, so the system (|13l) reduces to a 
system of 5 equations. 



C. Field equations 

The mean field values (/), a, luq and 69 are computed 
by solving their field equations, which are: 







4Box' lnx-2Bo5x' In 



Bo5x o b^Sx^ + e'lX 



2e\x rnlxiuol + (Aw^Aw^)) 

CTO 

mlx{bl + (A6^ • Ab^)) , 



(14) 



= 



.9^0 f 



[psp{-mli{(j)^, {A<j)f)), T) 



+ psu{nil,{^,,{A4>1)).p;.T)\ - [Boa^Sx ^ ^^—^ 



= AG\Lol + ml{x^ + {Ax^))ujQ-g^pB 
- m2(x' + (Ax'))5o".gpP3 



00 



v/A:2 -I- mf e''(V'^'+™f -K) _ 1 



where psp and pg^ are the scalar densities of protons 
(11) and neutrons, respectively, pj, = [pp — Pn)/2 and x — 
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FIG. 1: 3D-plot of the ratio cr/cro as a function of baryon density 
and temperature. Here Z/A=0.5 




FIG. 2: 3D-pIot of the ratio (j)/4>0 as a function of baryon density 
and temperature. Here Z/A=0.5 



wide range of temperatures, baryon densities and isospin 
asymmetries. One of the most interesting applications is 
the study of chiral and scale phase transitions. To this 
purpose we used the code to build a phase diagram, as 
an extension of Fig. 8 of Ref. pj, by computing (t/cto 
(see Fig. [1]) and ^/0o (see Fig. [2]) as functions of T and 
PB in the case of isospin symmetric matter. Notice that 
at large temperatures ct/ctq sharply decreases reaching a 
plateau where it assumes very small values. This plateau 
is a region where chiral symmetry is strongly restored, 
while scale invariance is still broken. When temperature 
becomes even larger scale invariance restores with a first 
order transition {(j>/4>o vanishes) [l^], dragging together 
the restoration of chiral symmetry {a/aQ vanishes). The 
code allows also to study the role of the isospin in chiral 
and scale phase transitions. 

It is necessary to point out that the results obtained 
here do not appear very different by those presented in 
Refs. (although some differences are observed in the 
isospin splitting of the pion mass). The thermal fluctu- 
ation of the dilaton field, indeed, does not play a funda- 
mental role. This is due to the value of the dilaton mass 
which remains large (about 1 GeV) also in the density- 
temperature region where the scale invariance restores. 

In the future it would be interesting to include the 
hyperons in calculations, through an extension of the 
model in SU(3). This direction is actually explored by 
the Frankfurt group (llj . The most interesting extension 
is the inclusion of the quark degrees of freedom, what 
will be done in the next works. 
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APPLICATIONS 



The CDM code [11,13] can be used for many applications References 
since it allows to compute many useful quantities in a 
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